Physical properties of noncentrosymmetric superconductor RU7B3 
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Transition metal boride RU7B3 was found to be a noncentrosymmetric superconductor with Tc equal to 3.3 
K. Superconducting and normal state properties of RU7B3 were determined by a self-consistent analysis through 
resistivityCp,-,- and p^) , specific heat, lower critical field measurement and electronic band structure calculation. 
It is found that RU7B3 belongs to an s-wave dominated single band superconductor with energy gap 0.5 meV and 
could be categorized into type II superconductor with weak electron-phonon coupling. Unusual 'kink' feature 
is clearly observed in field-broadening resistivity curves, suggesting the possible mixture of spin triplet induced 
by the lattice without inversion symmetry. 

PACS numbers: 74.70.Ad, 74.25.Qt, 74.25. Sv 



I. INTRODUCTION 

Unconventional superconductors have been extensively 
studied during the past decades for the underlying funda- 
mental physics or the potential industrial applications. Some 
well known examples are heavy fermion superconductors', 
cuprates^, Sr2Ru04^^ and the newly discovered iron arsenide'*. 
In recent years superconductors lacking of lattice inversion 
center have received intensive attentions for the possibility of 
spin-triplet dominated pairing symmetry. An important exam- 
ple is noncentrosymmetric superconductor CePt3Si^. Due to 
the nontrivial antisymmetric spin-orbit coupling(ASOC) ef- 
fect induced by the heavy Platinum atom and the lacking of 
inversion symmetry, unconventional superconducting proper- 
ties are observed, for instance the high upper critical field 
far beyond the Pauli-Clogston limit^. The subsequent nu- 
clear spin-lattice relaxation rate as well as magnetic penetra- 
tion depth measurements show line nodes in superconducting 
gap for CePt3Si^''^. For a material without inversion sym- 
metry the spin degeneracy is lifted by ASOC, under such 
a condition, orbital angular momentum(L) and spin angular 
momentum(5) are not good quantum numbers any more, thus 
the strict categorization of even-parity spin singlet and odd- 
parity spin triplet conformed to Pauli's exclusion and parity 
conservation is violated, then spin triplet may be mixed with 
spin singlet. Up to now, several noncentrosymmetric super- 
conductors have been reported, the development concerning 
such a scope could be consulted a overview given by Sigrist 
et. al.'. However strong electronic coiTelation in some mate- 
rials complicates the studies on ASOC effect, such as heavy 
fermion properties of CePtsSi. Recently the pairing symme- 
try of Li2Pt3B has been proved to be with line nodes'*^, such 
a material is not strongly correlated and could be regarded as 
an appropriate example for ASOC effect. For more compre- 
hension on ASOC, thus, any efforts searching for new non- 
centrosymmetric superconductors are worthwhile, especially 
in the absence of correlation effect. 

In the present paper we report a new noncentrosymmetric 
superconductor RU7B 3 with Tc equal to 3.3 K. Transport mea- 
surements as well as electronic band structure calculation gave 
a detailed description to the properties. As to our best knowl- 
edge, neither measurements nor calculations have been em- 



barked except Matthias mentioned only the Tc as 2.58 K in 
1950's". Additionally, for the first time we point out the non- 
centrosymmetric characterization of RU7B3. 

The paper is organized as follows: Section II describes 
the sample preparation process and the experimental details; 
then in section III we provide the structure illustration and 
the measurement of many transport properties, including AC 
magnetization, resistivity, magnetoresistance and Hall coeffi- 
cient; section IV focuses on the specific heat and lower crit- 
ical field measurements; in the following section V full po- 
tential electronic structure calculation is discussed; finally a 
self-consistent normal state and superconducting parameters 
are determined in section VI. 



II. EXPERIMENT 

Polycrystalline sample RU7B3 was synthesized by conven- 
tional solid state sintering process. Stoichiometric Ruthenium 
powder(purity 99.9%) and Boron powder(purity 99.95%) 
were mixed together and grounded thoroughly, then the mix- 
ture was pressed into pellet with a pressure 10 MPa. The pellet 
was wrapped with tantalum foil and sealed into an evacuated 
silicon tube. 95% high purity Ar mixed by 5% H2 was selected 
as protected atmosphere. The tube was warmed up to 1000 ° 
C in a muffle furnace with a rate 100 ° C/h and sintered for 
100 hours for composition homogeneity, then the furnace was 
cooled down to room temperature. The ingot was ceramic- 
like and with a silver gray color It was observed that the su- 
perconducting quality was stable and not sensitive to humid- 
ity or oxygen. X-ray diffraction (XRD) pattern measurement 
was performed at room temperature employing an M18AHF 
x-ray diffractometer (MAC Science). Cu Ka was used as the 
radiation target. Crystallographic orientation and index were 
determined by Powder-X'-, an software for processing x-ray 
diffraction data. The AC susceptibility were measured based 
on an Oxford cryogenic system (MagIab-Exa-12). The resis- 
tivity and specific heat were measured on the Quantum De- 
sign instrument PPMS with temperature down to 1.8 K and 
the PPMS based dilution refrigerator (DR) down to 50 mK. 
The temperatures of both systems have been well calibrated 
showing consistency with an error below 2% in the tempera- 
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TABLE I: Structure parameters of RU7B3 
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"No. represents the number of nearest Boron for different cites of Ruthe- 
nium. 

*/ stands for the interatomic distance between Boron ions and Ruthenium 
ions. 
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FIG. 1: (a) a-b projection of RU7B3 structure, the skeleton is con- 
sisted of metal tetrahedra and metal octahedra and interstitial Boron 
anions. From such a projection, the lattice with hexagonal rotation 
symmetry is obvious, (b) a-c projection of the lattice, two 'chains' 
are built up by ruthenium ions at different coordinations. 



ture range from 1.8 K to 10 K. The temperature dependence 
of magnetization was measured on the Quantum Design in- 
strument MPMS with temperatures down to 2 K. The lower 
critical field measurement was based on a two dimensional 
electron gas(2DEG) micro Hall probe with an active area of 
10x10 fj.np-. All M(H) curves were taken in zero field cooled 
mode with initial temperature up to 6 K. A low field sweep 
rate of 60 Oe/min was selected to measure isothermal magne- 
tization curves. Self-consistent band calculations were carried 
out using the linear muflintin orbital(LMTO) method on full 
potential plane wave representation(FP-LMTO-PLW)i^ 



III. STRUCTURE ILLUSTRATION AND TRANSPORT 
PROPERTIES 

The crystal structure of RU7B3 was determined by 
Aronsson'** in the late 1950's. It was found that the lattice 
is hexagonal with a space group P6j,mc. There are 20 atoms 
in one unit cell with eff'ective coordinate Ru/(6c), Ru//(6c), 
Ru///(2b) and B(6c), respectively. Thus two formula units ex- 
ist in one unit cell. In that paper the author gave a relative 
comprehensive description to the lattice structure of RU7B3, 
however, thanks to the complex structure of transition metal 
boride, a more detailed illustration should be added as a sup- 
plement for better comprehension on it's properties. As to 



space group P6}mc, a distinct characterization of the crystal 
lattice is without inversion symmetry, for example, the cen- 
ter of Boron atoms' sub-lattice dislocates the counterpart of 
Ruthenium atoms along c-axis, thus the inversion symmetry 
is broken along such direction. For more clear understand- 
ing on the crystal lattice, we illustrate the structure along two 
types projection in Fig. 1. Fig. 1(a) shows the ab projec- 
tion of RU7B3 structure. It is found that the skeleton consists 
of metal tetrahedra and metal octahedra and interstitial Boron 
anions. From such a projection, the lattice with hexagonal 
rotation symmetry is obvious. Fig. 1(b) is the ac projec- 
tion illustration. A very interesting phenomenon is that two 
'chains' are built up by ruthenium ions at different coordina- 
tions. The metal octahedra at each corner of the lattice is built 
up by Ru//(6c), those octahedras share face along c-axis and 
thus a zig-zag chain is formed as shown in Fig. 1(b). The 
rest Ru/(6c) and Ru///(2b) form two tetrahedras in one unit 
cell at different (x,y) positions, along c direction two tetrahe- 
dras (at the same (x,y) positions) share face and then forming 
a hexahedra, each hexahedra is connected by Ru///(2b) and 
thus another type of row is formed. The special structure con- 
figuration might play an important role in transport properties, 
also the environment of Ruthenium ions (including the inter- 
atomic distance and the number of the nearest boron ions) is 
very important. Detailed parameters are included in Table 1 
as shown below. 

Fig. 2 shows the XRD pattern of the sample RU7B3, which 
can be indexed in a hexagonal symmetry with a-b - 7.4629 
A and c - 4.7141 A. The indexed indices slightly deviate from 
the reported parameters a - b - 7.467 A and c - 4.713 A. It 
is clearly found that most diffraction peaks are indexed except 
one minor peak possibly from un-reacted boron. From the 
quality of XRD data, it is estimated that the purity of RU7B3 
we synthesized is about 95%. In the following specific heat 
measurement, 90% superconducting component also prove 
the sample's purity. 

Fig. 3 presents the AC susceptibility measurement under 
different magnetic fields. It is found that a sharp diamag- 
netic transition happens at 3.3 K and the diamagnetic value 
approaches a constant at 3K, 0.3 K transition width indicates 
the good superconducting quality. Applying magnetic fields, 
the transition curve moves parallel to low temperatures. For 
estimating the upper critical field, we take 95% diamagnetic 
value as criterion and will discuss the result in the following 
paragraph. 
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FIG. 2: (Color online) X-ray difFraction of RU7B3, most diffraction 
peaks are indexed except one minor peak possibly from un-reacted 
boron, the purity is estimated about 95%. The asterisk marks the 
peak from impurity phase. 

Fig. 4 shows the resistivity data from 2 K to room temper- 
atures, the curve shows a good metallic behavior with a zero 
temperature resistivity(po) 9 yuQ-cm. Such a high conductivity 
is the common feature of transition metal boride. It is found 
that the residual resistivity ratio(RRR) is 28, the relative high 
RRR in polycrystalline sample indicates that impurity scatter- 
ing is trivial for conductivity. The inset of Fig. 4 shows the en- 
largement of superconducting transition, the resistivity drops 
sharply to zero at 3.3 K. Thus the resistivity data and AC sus- 
ceptibility give a self-consistent Tc=3.3 Kfor RU7B3, which is 
slightly higher than the preliminary mentioned rc~2.58 K by 
Matthias" . Considering the high purity of the sample, a quan- 
titative analysis of normal state resistivity is deserved. We try 
Wilson's model for transition metals'^ 



FIG. 3: (Color online) AC susceptibility of RU7B3, a sharp super- 
conducting transition happens at 3.3 K under zero field, a narrow 
transition width less than 0.3 K indicates the good superconducting 
quality. 
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where is a cutoff similar to the Debye temperature and Ksd 
is a constant. Using parameters ©5^=500 K and Ksd-1350, a 
good fitting is obtained as shown in Fig. 4. 

Fig. 5 shows the field-broadening resistivity measurement 
on polycrystalline RU7B3 down to 100 mK. When magnetic 
field is applied, the onset parts of transition is rounded very 
similar to the feature of critical fluctuation in superconduc- 
tors. While adding fields up to 0.6 Tesla, the rounded part 
of transition evolves into a 'kink' structure as shown clearly 
in Fig. 5. It seems that the kinks break the superconduct- 
ing transition curves into two parts, the lower part moves 
quickly to low temperatures showing a strong dependence of 
magnetic fields, while the upper parts respond reversely. It 
is observed that the zero resistivity point approach zero tem- 
perature when sample is bearing LIT magnetic field, while 
superconductivity is depressed completely at about 5 T. It is 
thus very interesting to note that the determined upper critical 
field strongly depends on the selected resistive criterion. As 



FIG. 4: (Color online) Resistivity of RU7B3 under zero field is mea- 
sured from 2 K to room temperatures, a good fitting is given by Wil- 
son s-d scattering model. The inset shows the enlargement of super- 
conducting transition, the transition width is less than 0.3 K. 



shown in Fig. 6, the criterion of 99% p„ and zero resistivity 
will lead to two distinct Hc2(0) with a ratio of about 5. A 
conventional interpretation, as discussed in MgioIrigBe"', for 
field broadening transition is filamentary-like superconductiv- 
ity along grain boundaries, the stronger scattering reduces the 
mean free path and consequently influences on the coherence 
length. However, we argue that unlike the case of MgioIrigBe 
the clear kink feature in RU7B3 complicates the determination 
for Hc2, also the kink doesn't come from sample inhomogene- 
ity or two phases because of the sharp transition in low fields. 
As to our knowledge, the unconventional 'kink' or 'step' fea- 
tures in transition curves have been extensively discussed in 
single crystals MgB2. Several reasons have been attributed 
to this eff'ect, including superconducting fluctuation'^, two 
superconducting gaps'^, surface barriers'"^ and vortex lattice 
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FIG. 5: (Color online)Field broadening superconducting transition 
curves of RU7B3 from 0. 1 K to 4 K, a clear kink feature appears as 
the field exceeds 0.6 T. The drop of resistivity is totally suppressed 
by 5 T magnetic field. Distinct magnetroresistance is observed. 

meltin g^°i^' ,etc. However taking account of the very low Tc 
and Hc2, superconducting fluctuation is believed to be weak. 
For surface barriers and vortex lattice melting, the polycrys- 
talline quality seems to exclude their possibilities. As to the 
two-gap scenario, the following specific heat and lower crit- 
ical field measurements oppose such a point of view. Thus 
we assumed the possibility that spin triplet induced kink fea- 
ture in the framework of inversion symmetry is broken. It is 
known that applying fields broke time reversal symmetry and 
is detrimental to spin singlet, whereas triplet pairing remains 
unaff'ected. So for a superconductor with pairing symmetry 
mixed by singlet and triplet, the kink feature in field broad- 
ening resistivity curves could be plausible on certain extent. 
In Fig. 6 we plot the phase diagram of RU7B3, the criteria 
are taken as below, for resistivity 99%p„ and zero resistance, 
for magnetization 95% diamagnetic signal and the half po- 
sition of specific heat anomaly for thermodynamic measure- 
ment. Thus the derived AHcildT equal to -0.43 T/K or -0.277 
T/K for criteria 99%p„ and zero resistance, respectively. It is 
found that except 99%p„ the other three criteria determined 
data points overlap almost together Further consideration is 
that bulk properties provided by specific heat measurement, 
thus in the following discussion we use 1.1 T as Hc2(0) and 
-0.277 T/K as dHci/dT near Tc. We also try to use Ginzburg- 
Landau formula to fit the data points determined by 99%p„, 

HaiT) = //c2(0)— 2 , (2) 
1 + f 

where t is the normalized temperature T/Tc, it is found that 
the fitting curve strongly deviates the 99%p„ points as shown 
in Fig. 6, indicating the invalidity of Landau twice order phase 
transition theory in the present material. 

Another distinct characterization of Fig. 5 is the field 
induced magnetoresistance. Typically magnetoresistance is 
used to investigate the electronic scattering process and pro- 



FIG. 6: (Color online) Phase diagram of RU7B3, a non-trivial flux 
flow area for criteria 99%p„ and zero resistance is formed by kink 
shape of transition curves. 95% diamagnetic signal from AC suscep- 
tibility and one half points of specific heat jumps are also selected 
as criteria to estimate the intrinsic value of HciiO)- It is shown that 
except for the criterion of 99%p„ the other points from different cri- 
teria overlap together, thus the experimental data 1. 1 T is determined 
as Hc2(0)- It is found that the estimation of Hc2 from WHH formula 
fits experimental data very well, whereas Ginzburg-Landau formula 
fails to describe the upper bound of Hc2- 



vide useful information on fermi surface (FS) topology. So 
detailed studies are needed. In Fig. 7 we present the tem- 
perature dependence of resistivity under magnetic fields from 
to 9 T. It is found that the magnetoresistance(MR) is about 
16% [(pgr - PotVpot] 3.1 5 K, which is one order of magni- 
tude larger than the ratio of recently discovered iron arsenide 
LaFeAsOi-.vFj^^. The latter was regraded as superconductor 
with multiple bands. A simple verification for the possibility 
of multigap effect is the scaling based on Kohler's rule. The 
Kohler's rule is written as Ap/po=F(H/po), where F is a func- 
tion depending on the nature of the metal itself. For single 
band metal with symmetric Fermi surface topology Kohler's 
law should be conserved. It is shown as Fig. 7(b) that Kohler's 
rule is only slightly violated. Unlike typical multi-band su- 
perconductor MgBi^^ and LaFeAsOi-vF.i, the breakdown of 
Kohler's law is trivial in RU7B3, indicating a dominated sin- 
gle band behavior. The further specific heat and lower crit- 
ical field measurements further provide the same conclusion. 
However we believe that the slightly violation of Kohler's rule 
could be induced by noncentrosymmetric structure of RU7B3, 
due to ASOC, the degenerate spin-up and spin-down bands 
are spilt, so Kohler's rule is slightly broken. 

Hall coefficient (Rh) measurement was done by sweeping 
temperature at magnetic field 9 T and reversing field (-9 T). 
For avoiding the possible temperature hysteresis, increasing 
temperature mode with a moderate rate 1 K/min was adopted 
for both positive and negative fields. The Hall coefficient is 
shown in Fig. 8, it is found that the charge carrier of RU7B3 
is dominated by hole-like carriers with R// 3~6xlO"''^ m^/C 
from 2 K to 200 K. For verifying the Rh, we also use the val- 
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FIG. 8: (Color online)Hall coefficient measurement of RU7B3 from 
sweeping temperature under reversing fields, the charge carrier is 
dominated by hole from low temperature to 200 K, the density of 
charge carrier is about 1 x 10^^/cm"^. The squares represent the data 
measured by sweeping fields at fixed temperatures. 



IV. SPECIFIC HEAT AND LOWER CRITICAL FIELD 
MEASUREMENT 



FIG. 7: (Color online) (a) Temperature dependence of resistivity of 
RU7B3 under magnetic fields from to 9 T, (b) the derived mag- 
netroresistance Ap/po is 16% at 5 K, it is shown that the scaling of 
Kohler's law is slightly violated. The violation is attributed to band 
splitting by ASOC instead of multiband effect. 



ues derived from sweeping field at three temperature points 
2 K, 100 K and 200 K. The low temperature Rh is consis- 
tent with the value from sweeping temperature, while error 
bars exist in high temperatures. The charge caiTier density 
calculated by n=l/(R//-e) is about IxlQ^^ / cm^ , which is two 
order of magnitude larger than low carrier density supercon- 
ductors, for example, cuprates^'* and hole doped iron arsenide 
(Lai-j^Sr JFeAsO^^. We also notice the nonlinear tempera- 
ture dependence of Rh, however, the relative change of Rh 
from 2 K to 200 K is small as shown in Fig. 8. Moreover, 
it is known that the hall effect is very sensitive to the tem- 
perature dependent scattering rate, local fermi velocity and 
complex FS topology^^, thus considering the polycrystalline 
quality detailed analysis on temperature dependent Hall coef- 
ficient would not be discussed here. 



A. Specific heat 

Fig. 9 shows the raw data of specific heat under differ- 
ent magnetic fields from zero to 3 T. With increasing fields 
the specific heat anomaly near Tc move quickly to low tem- 
peratures leaving a background consistent with that above 
at zero field. Thus the normal state specific heat could be 
extracted easily with the relation C/T - y„ + ySr^, where 
y„ is the normal state specific heat coefficient and /? coiTe- 
sponds to phonon contribution. It is found that y6=0.3732 
mJ/molK'^ and y„= 89.95 mJ/molK', however a residual 
value yo ~ 9.8 mJ/molK^ indicates the existence of about 
10% non-superconducting fraction. The non-superconducting 
fraction could partly come from unreacted boron as inferred 
from analysis from X-ray diffraction. Thus the normal state 
Sommerfeld constant could be determined from the rela- 
tion je-jn-Ji'i as 80.15 mJjmolK}. Using the relation ©£> 
=(12;r'*;tBA^AZ/5y6)'^^ where A^a = 6.02x10^3 the Avogadro 
constant, Z=20 the number of atoms in one unit cell, we get 
the Debye temperature ©o = 470.18 K. In the previous sec- 
tion, we obtained a similar value 0sj=5OO K from the resis- 
tivity curve fitting, the consistent values prove the reliability 
of two different measurements. It is noticed that comparing 
with the values of MgiolrigBg^^ and Li2Pt3B^^ 7„ in our mea- 
surement is relatively high, thus a prudent checking is neces- 
sary. For a type-11 superconductor, y„ could be estimated as 
the following relation^^: 
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FIG. 9: (Color online) Raw data of specific heat plotted as C/T vs. 
T^. All filled symbols represent the data taken with the DR based 
on the PPMS at various magnetic fields. The open squares show 
the data taken with the PPMS at different fields. The thick solid 
line represents the normal state specific heat which contains both the 
phonon y,,i, and the electronic contributions. 



where A = 3me/n^kB^4-479.2HT /K)(nm)-\j/m^K^)-\ 
using -dnoHc2(T)/dT\T, « 0.277 T/K, p„ = 9 yuQcm, and 
taking t] - I for the weak coupling case, we have y„ = 94 
mJ/molK^ which is very close to the upper bound of the ex- 
perimental value 89.95 mJ/molK^. Further we could estimate 
the Ae-ph via McMillan equation-': 



Tr = 



1.45 



exp(- 



1.04(1 +^,_P,) 



•|/*(1 +0.62^_p^) 



(4) 



where fi* is the Coulomb pseudopotential taking 0.11, 
©0=470.175 K and Tc=3.3 K, we obtain ^,_p,,=0.48. The 
value indicates that RU7B3 belongs to a weak coupling super- 
conductors. 

For noncentrosymmetric superconductors, novel pairing 
symmetry could be achieved due to the mixing of spin singlet 
and spin triplet. Specific heat is a useful tool to investigate 
the material's low energy excitation. So, subtracting the con- 
tribution of phonon, we present temperature dependence of je 
under magnetic fields up to 3 T in Fig. 10(a). For the conve- 
nience of theoretical analysis, we further subtract y„ of zero 
field data as shown in Fig. 10(b). Thus the weak coupling 
BCS formula could be used: 
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where f - + A^{T, 6). In obtaining the theoretical fit 
we take the implicit relation Ao(r) derived from the weak 
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FIG. 10: (Color online) (a) Temperature dependence of for mag- 
netic field up to 3 T, (b) Temperature dependence of y„e-ye at zero 
field. The blue solid, green dashed and orange dotted lines are theo- 
retical curves calculated based on BCS model with a gap of s wave, 
d wave and p wave, respectively. 



coupling BCS theory for superconductors with different pair- 
ing symmetries: A(r, ff) - A()(r) for s-wave, A{T, 9) - 
Ao(r)cos20 for d-wave, and A(r,0) = Ao(r)cos0 for p- 
wave, respectively. The theoretical curve of s wave fits the 
experimental data very well leading to an isotropic gap value 
Ao = 0.5 meV and = 3.22 K. The ratio Ao/keT, = 1.80 ob- 
tained here is quite close to the prediction for the weak cou- 
pling limit (Aq/ ksTc - 1.76). In addition, the specific heat 
anomaly at T^. is ACe/y„T\T^. x 1.31 being very close to the 
theoretical value 1 .43 predicted for the case of weak coupling. 

Condensation energy(£'c) is an important parameter for su- 
perconductor, thus we calculate Ec with the following pro- 
cess, firstly the entropy difference between normal state and 
superconducting state could be obtained by 5 „ - 5 , = (y„ - 

ye)dT' , and then Ec is calculated through Ec - j^'^{S„ - 
S s)dT'. The resulted temperature dependence of Ec is shown 
in Fig. 11, the inset is the entropy difference between normal 
state and superconducting state. The Ec is about 192 mJ/mol 
at K. Alternatively, Ec could be calculated by the following 
equation: 



Ec = aN(EF)Al/2 = a 



3 1 
4^^ 



(6) 



For a BCS s-wave superconductor, a^l, taking y„=80.15 
mJ/mol-K^ and Ao=0.5 meV, we obtain a value 205 mJ/moI, 
which is close to experimental value 192 mJ/mol. In addition. 
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FIG. 11: (Color online) Superconducting condensation energy of 
RU7B3 calculated by specific heat, the inset shows entropy difference 
between the normal and superconducting state. 



the consistence reversely verifies the validity of 7„ and Aq de- 
termined through our experiment. From condensation energy, 
the thermodynamic critical field yuoHc(O) could be calculated 
via the relation /io//^(0)/2 - Fm-Fs - ffijn - 7e)dT, yield- 
ing jUoHc(0)= 612 Oe. For a comparison, //oHc(O) of another 
noncentrosymmetric superconductor MgioIrigBe is about 300 
Oe'^ 



FIG. 12: (Color online) The raw data of M(H) curves at different 
temperature. The red dashed line give the Meissner linear approach. 
It is found that the value of M is very approaching to that of mag- 
netic field, indicating almost 100% diamagnetization at lower tem- 
peratures. 



bar due to noise induced uncertainty. According to BCS the- 
ory for clean superconductors, the normalized Hci(T)/Hci(0) 
is expressed as follows^: 



B. Lower critical field measurement 

Lower critical field(Hci ) is an important parameter for a su- 
perconductor According to the Ginzburg-Landau theory, Hci 
reflects the superfluid density ps since Hci is related to Lon- 
don penetration depth A and thus a relation is established that 
Hci~l//l^. Moreover, the temperature dependence of Hci, es- 
pecially the low temperature features, is always used to inves- 
tigate the superconducting pairing symmetry and multi-gap 
effect, for instance the node feature for pairing symmetry and 
gaps' value for the latter In this section we used a two di- 
mensional electron gas(2DEG) micro Hall probe to measure 
local magnetization loops of RU7B3. For weakening the com- 
plex effects of the character of field penetration, such as Bean- 
Livingston surface barriers and geometrical barriers, we used 
a low field sweep rate of 60 Oe/min to measure isothermal 
magnetization curves. 

Fig. 12 is the initial isothermal M(H) curves over the tem- 
perature range from 1.22 K to 3.2 K. It is found that the low- 
field parts of those M(H) curves overlap almost on one line 
(red dash line is guided to the eyes in Fig. 12) with a con- 
stant slope, which is attributed to Meissner effect and called 
as Meissner line. Thus Hci could be determined as the de- 
parture between M(H) curve and Meissner line with the same 
criterion for all curves. The resulted temperature dependence 
of Hci (T)(normalized to Hci(O)) is shown in Fig. 13. The 
inset of Fig. 13 shows the criterion for determination of the 
value of Hci at 1.4 K, a value of 10 Oe was regarded as error 



Hem m 



= 1-2 (- 

Ja(T) 



df{E) 
d{E) 



)D{E)d{E), (7) 



Where A(r) is the BCS superconducting energy gap, 
f{E)=l/[e\p{-E/kBT)+l] is the Fermi distribution function, 
and D(E)=E/[E^-A^(T)]'^^ is the quasiparticle density of 
states. We use above equation to fit the experimental data, 
Hci(O) and A(0) are fitting parameters. It is found that single 
gap s-wave pairing could give an appropriate fit with fitting 
values Hci (0)= 110 Oe and A(0)=0.5 meV, the latter is con- 
sistent with that of specific heat measurement. Thus the good 
consistence indicates the reliability of results determined by 
Hci measurement, although the lower temperature (less than 
1.2 K) experimental data is absent. A possible argument is 
that the nominal Hci obtained from experimental data fitting 
may not reflect the true value due to the inevitable surface 
baiTier and geometrical baiTier induced by the polycrystalline 
quality. Thus in the following paragraph we will do self- 
consistent checking from superconducting parameters deter- 
mined by other measurements. 



V. ELECTRONIC STRUCTURE CALCULATION 

In this section we present density of states(DOS) and band 
dispersion results based on full potential linear-muffin-tin- 
orbital program LMTART by Savrasovl^;^. Full potential 
approximation. Plane wave expansion(PLW), is selected and 
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FIG. 13: (Color online) The extracted Hc\ as a function of temper- 
ature, isotropic s wave with gap value 0.5 meV could give a good 
fitting to the experimental data. The inset shows the criterion for Hc\ 
determination at temperature 1 .4 K. 
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FIG. 14: (Color online) Density of states of RU7B3 calculated by full 
potential PLW-LMTO, the Dos at chemical potential is about 20.988 
state/(eV cell spin). It is shown that electronic structure is dominated 
by Ruthenium 4d states, boron 2p state contributes weakly. 



believed to give the adequate accuracy. Fig. 14 is the DOS 
calculation of RU7B3. The total DOS curve have numerous 
van Hove singularities, the feature is very similar to that of 
MgioIrigBie^^, in that paper the author attributed the charac- 
terization as large numbers of atom in the unit cell and var- 
ious interatomic distances. In RU7B3 only 20 atoms exist in 
one unit cell, thus the numerous van Hove singularities could 
stem from the various interatomic distances complicated by 
the lacking of a inversion center. Another feature is that elec- 
tronic structure is dominated by Ruthenium 4d states, boron 
2p state contributes weakly. It is reasonable to understand 
from structure aspect that the lattice of RU7B3 is mainly con- 
sists of metal tetrahedra and metal octahedra or 'chains' along 
c-direction, thus charge carries naturally favor those special 
channels in a crystal lattice. The total DOS at chemical poten- 
tial for RU7B3 is 20.988 state/eV per formula unit. For check- 
ing the calculated DOS at chemical potential, we could simply 
estimate the DOS from jn in the framework of free electron 
gas. 
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FIG. 15: (Color online) Band dispersion curves of noncentrosym- 
metric materials RU7B3, the bands are doubly bound due to the 
ASOC induced spin-up and spin-down splitting. 



A^(0) = (^r'-r«, (8) 

where j„ is selected as 80.15 mJ jmolK^ , k^ is Boltzmann 
constant=l. 380658x10"^^ J/K, N(0) represents the density of 
state. The obtained N(0) is about 17 state/cell per formula, 
which is close the calculated value 20.988. Fig. 15 is the 
band dispersion curves near Fermi energy. A distinct feature 
is that the all bands are doubly accompanied, the feature is 
due to asymmetric spin-orbit coupling effect, thus the degen- 
eracy of spin-up and spin-down is lifted. It is noticed that 
at some k-points with high symmetry splitting instead of de- 
generacy also exists, which could come from problems such 
as un-adequate optimized parameters at initialization during 
computation. Nevertheless, it is believed that such stigma can 
not affect the main results of the present paper. 



VI. SELF-CONSISTENCE AMONG SUPERCONDUCTING 
PARAMETERS 

In the discussion of magnetoresistance, the Kohler's rule 
was found to be only slightly broken, indicating a symmetric 
Fermi surface topology. Thus we could deduce the Fermi- 
wave number(A:/r) from charge carrier density(n)-''* assuming 
a single spherical Fermi surface, fef =(3;7r^n)'^^=6.6527 nm"', 
where n=0. 995x10^^ cm"^. The effective mass is estimated 
as m* - i3Ti^y„)/(V,„oik^gkf)- ITm^,/, where m^,/ is bare elec- 
tron mass and molar volume y,„o/=136.89 cm^/mol. Then the 
Fermi velocity Vf-h kf/m* is about 0.47x10^ m/s. The mean- 
free-path is evaluated as l-fikf I {pane^)-^ 1 -36 nm. The super- 
conducting penetration depth A(Q) - ^jr^F]Ji^ne^ is 214 nm. 
So the coherence length could be estimated using the BCS ex- 
pression ^(0) = Q.\%tiVfl{kBTc) as 19.5 nm-". Thus the above 
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TABLE II: Superconducting and normal state properties of RU7B3 



parameters 


RU7B3 


7„(mJ/mol K-) 


80.15 


A'(0)(state/eV cell spin) 


20.988 


A (m eV) 


0.5 


Ha (Oe) 


11000 


Hci (Oe) 


90 


f (nm) 


17.3 


A (nm) 


215 


K 


12.4 


HciO) (Oe) 


612-628 


A/kfiTc 


1.80 


A,/r„Tc 


1.31 


/3 (mJ/mol K") 


0.3735 


0D (K) 


470.175 


£f (mJ/mol) 


192-205 




17mp 


n (cm"') 


1x10^2 


/(nm) 


31.36 


p„ (;^n- cm)(4 K) 


9 



superconducting parameters could give a stringent checking 
on experimental data, such as ^(0) and Hci and ^oHciQ). In 
the upper critical fields measurement, Hc2(0) is determined 
as 1.1 T, so using ^(0)= ^J(f>Q/i2nHc2{0)), where 0o is flux 
quanta, coherence length is 17.3 nm, such a value is very close 
to the deduced ^(0) 19.5 nm. For checking on Hci , the follow- 
ing formula is used'^; 




60 90 

T(K) 



150 



150 



FIG. 16: (Color online) Normal state magnetization measurement of 
RU7B3 from 5 K to 150 K. (a) FC(field cooling) mode, a rounded 
hump happens at about 40 K, (b) ZFC(zero field cooling) mode, the 
rounded hump on FC curve changes into peak and moves to higher 
temperature about 55 K. 



(9) 



yielding fioHci = 90.2 Oe, where Aq and ^0 are deduced values 
from charge carrier density. In our lower critical field mea- 
surement we obtain jUoHci = 110 Oe, which is larger than the 
estimated value 90.222 Oe. The thermodynamic critical field 
jUoHc(O) could be obtained from the following formulai^: 



HciHc2 = Hl(0)lnA 



(10) 



using Hc2=l 1000 Oe, Hci =90.2 Oe, /lo=214 nm and ^o=17.3 
nm, /voHc(O) is given as 628 Oe, which is very close to 
the value of 612 Oe determined by specific heat. For fur- 
ther checking on the experimental value of Hci (0), firstly the 
experimental values //ci(110 Oe) and ^(0)(17.306 nm) are 
used into Equation 9, the solved /lo= 189.3 nm, then using 
that value Hc(0) could be estimated by Equation 10, yield- 
ing jL/oHc(0)=71 1 Oe, which is about 100 Oe larger than that 
obtained from specific heat measurement. Therefore it is safe 
to conclude that the intrinsic value of Hci(O) is about 90 Oe. 
So in Table II we list the superconducting and normal state 
properties of noncentrosymmetric material RU7B3. 

In the final part of self-consistent checking on physical pa- 
rameters of RU7B3, a brief discussion on the possible exotic 



properties due to ASOC is necessary. Generally, two crite- 
ria for novel superconductivity in the framework of noncen- 
trosymmetry have been established. The first one is Pauli- 
Clogston limiting field*", which could be expressed as Hp(Q) = 

A(0)/2V2yUB = 1.83rc, for RU7B3 with Tc 3.3 K, yield- 
ing Hp{Q) X 6 T. In our experiment the upper bound of 
Hc2(99%p„) is found to be about 5 T, the experimental value 
is obvious below Pauli limit, indicating that Hc2 is still de- 
termined by orbital depairing fields. Nevertheless, the ob- 
served kink feature is still a puzzle under the common think- 
ing, such as two bands or superconducting fluctuation scenar- 
ios, et. al. Another criterion is the presence of line or point 
node in the superconducting gap. The good fitting of specific 
heat data with isotropic s-wave has provided a strong evidence 
that spin singlet dominated the condensate. The lower critical 
field measurement also give the same conclusion even though 
the lower temperature data (less than 1.2 K) is unfortunately 
absent. Thus a safe conclusion could be given that noncen- 
trosymmetric superconductor RU7B3 is dominated by s-wave 
pairing symmetry, minor spin triplet could admix among the 
majority of spin singlet. 

In the end of experiment and discussion, we show magne- 
tization of RU7B3 from 5 K to 150 K bearing the magnetic 
field modes of FC(field cooling) and ZFC(zero field cooling) 
in Fig. 16 . It is found that there is a hump at 40 K on the FC 
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curves, while for ZFC mode the hump changes into a peak and 
moves to higher temperature at about 55 K. For now it is still 
hard for us to comprehend the normal state magnetization of 
RU7B3, especially more subtle features adding on ZFC curves 
as shown in Fig. 16(b). For comparison we measured room 
temperature magnetization of Ruthenium element(not shown 
here), it is found that Ruthenium element is paramagnetic at 
300 K, lowering temperature an antiferromagnetic transition 
happens at 150 K on the paramagnetic background. The den- 
sity states calculation have shown that ruthenium contributes 
most DOS on FS. Thus from such a point view, the normal 
state magnetization of RU7B3 could be similar to that of ruthe- 
nium element. Furthermore for the ruthenium element the 
outer shell electron configuration is 4(f5s^, thus in RU7B3 the 
ruthenium cations with higher spin angular momentum could 
be anticipated. However, we could not exclude the possibil- 
ity of impurity induced magnetism due to the 10% residual 
yo for specific heat measurement. Nevertheless, the studies 
on magnetism of RU7B3 is worthwhile, for example the pos- 
sible antiferromagnetic fluctuation induced superconductivity 
has been a hot issue in MgCNi3^^ for RuSr2GdCu20g^*' the 
interplay between ferromagnetism and superconductor is also 
very attractive. 

VII. SUMMARY 

Transition metal boride RU7B3 was found to be a noncen- 
trosymmetric superconductor with Tc equal to 3.3 K. Super- 



conducting and normal state properties of RU7B3 were de- 
termined by a self-consistent analysis through the results of 
resistivity(p^;c and Pxy) , specific heat, lower critical field mea- 
surement and electronic structure calculation. It is found that 
RU7B3 belongs to an s-wave dominated single band super- 
conductor with energy gap 0.5 meV and could be categorized 
into type II superconductor with weak electron-phonon cou- 
pling. Unusual 'kink' features are clearly observed in field- 
broadening resistivity curves, possibly indicating the admix- 
ture of spin singlet and spin triplet due to the absence of lattice 
inversion symmetry. 
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